Dilepton excess from local parity breaking in baryon matter 
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Abstract 

^ cj We propose that local parity breaking induced by a large-scale fluctuation of topological 

O-i charge at large temperatures and/or condensation of pseudoscalar mesons in the isotriplet chan- 

nel for large baryon densities may be responsible for the substantial dilepton excess that is found 
for low invariant masses and moderate values of pt in central heavy ions collisions. This insofar 
■ unexplained enhancement could be understood by a combination of two effects leading both 

to an excess of e + e~ and pairs with respect to theoretical predictions based on conven- 

er - ) ■ tional hadronic processes: (a) a modification of the dispersion relation of photons and vector 

^ ' mesons propagating in such a nuclear medium due to local parity breaking; (b) the appearance 

m. of new decay channels, forbidden by parity conservation in QCD in the usual vacuum. Possible 

signatures of this effect and perspectives for its detection are discussed. 
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1 Introduction 



The possibility of local parity breaking (LPB) in baryon matter at extreme conditions has been 
discussed recently [IJEIEIII] in association with relativistic heavy ion collisions (HIC) at RHIC, GSI, 
and CERN [5]. It has been suggested in pQ that at high temperatures an isosinglet pseudoscalar 
background could arise due to large-scale topological charge fluctuations (studied recently in lattice 
QCD simulations [6j). These considerations led eventually to the formulation of the so-called Chiral 
Magnetic Effect (CME)[1J, claimed to have already been detected in the STAR and PHENIX 
experiments at RHIC [TJ. The effect should be most visible for noncentral HIC where large angular 
momenta induce large magnetic fields contributing to the chiral charge separation. In a separate 
development it was proposed in [3] that the presence of a phase where parity was spontaneously 
broken could be a rather generic feature of QCD at finite density. This conclusion was reached 
using effective Lagrangians; a proof of principle is unfortunately still missing in QCD due to the 
difficulties of dealing with non- vanishing chemical potentials in lattice QCD. 

On the other hand, in central HIC an abnormal dilepton excess has been observed [H [9] (past 
experiments are reviewed in [5]). In the PHENIX experiment, for instance, an excess is seen for 
dileptons in the range M < 1.2 GeV, for centrality + 20% and px < 1 GeV [9]. Theorists have so 
far been unable to account for this dilepton excess convincingly. Thermal effects inducing vector 
resonance broadening and/or mass dropping fall short of providing a full explanation [10j-|14j. 
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The dilepton excess is conspicuously absent for peripheral HIC (where the CME should be more 
visible). We conjecture that the CME and the 'anomalous' dilepton excess may be complementary 
effects revealing two facets associated to the formation of a thermodynamic phase where parity is 
locally broken. This is the main contention of this work. 

We shall assume the emergence a of time-dependent pseudoscalar background (associated with 
axial chemical potential induced by a non-trivial topological charge [TJ [2] or with pseudoscalar con- 
densation |3|) in central HIC0 and search for possible manifestation of local parity breaking (LPB) 
in electromagnetic probes. We investigate whether these changes may induce a large dilepton ex- 
cess in certain HIC. Without excluding the relevance of more conventional mechanisms to partially 
account for the enhancement, we shall conclude that a good part of the excess of dileptons might 
be a consequence of LPB due to the generation of a pseudoscalar, isosinglet or neutral isotriplet, 
condensate whose precise magnitude and time variation depends on the dynamics of the HIC. 

This paper is mostly concerned with the possibility of LPB associated to the appearance of 
local fluctuations of the topological charge (or equivalently, as we shall argue below, of a axial 
chemical potential). The issue of isotriplet condensation induced by finite densitjH will not be 
discussed here (see [3] for some tentative considerations). As we will see, this form of LPB can be 
incorporated by adding a parity-odd term in the Lagrangian; however, this breaking, while being 
technically 'soft', can be numerically large, leading to important consequences. We shall review in 
detail how in-medium hadronic resonances are influenced by this effect and discuss possible ways 
to verify its existence. We will briefly discuss how hadronic processes are substantially modified 
by the presence of the parity breaking piece and point out their likely relevance for the dilepton 
enhancement. However a detailed discussion of the latter aspects as well as a complete comparison 
with the experimental data are postponed to a subsequent publication. 



2 Generation of axial baryon charge and axial chemical potential 

Let us outline the relationship between emergent topological charge, baryon axial charge and axial 
chemical potential. Let us assume that a jump of topological charge T5 localized in a finite volume 
arises (for instance, in a hot environment due to sphaleron transitions [161 117j ) as a consequence of 
a HIC, and survives for a sizeable lifetime in the fireball, At ~ r fireball — 5 — 10 fm. For a given 
gauge field G^(x) the topological charge itself 

T sW = ^2 J 1 d * xK ^ K n = ^W- Tr (cFdPG* - i^G u GPG^ , (1) 

is not gauge invariant under large gauge transformations. Conventional notations for the gluon 
field and the gluon field-strength are used [I]. However, the jump in T5 can be associated to the 
space-time integral of the gauge-invariant Chern-Pontryagin density 

at 5 - m f) - r, ( o) - JL, f' dt j^ Atv ( G -g„„) = JL f' dt j^ * xlf . Kr . (2) 

We shall neglect here the topological current flux through the fireball boundary during the ther- 
modynamic phase. 

We could refer the reader to [15] where it has been shown that a pseudoscalar field slowly evolving in time dras- 
tically changes the electromagnetic properties of the vacuum. In particular a photon propagating in this background 
with sufficiently high energy may decay on shell in medium into dileptons. 
2 Finite density may contribute to isosinglet pseudoscalar condensation too. 
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It is well known that for the color- and iso-singlet axial current J^^ = q^^^q the local PCAC 
relation is afflicted with a gluon anomaly 

N f 

d 1 "^ - 2iqm qlb q = l^^K^. (3) 

This exact equation allows to associate a non-zero topological charge with a non-trivial axial charge 
Q\. Indeed, one can integrate the local PCAC ([3]) over a finite space volume of fireball where a 
nontrivial topological charge is located and find the connection between their time derivatives 

-(Q q 5 -2N f T 5 )~2i [ d 3 xqm gl5 q, Q\=l d 3 xq lQlb q. (4) 
at ./vol. ./vol. 

In the previous equation the fluxes across the fireball surface of the axial and topological currents in 
the absence of strong magnetic fields are neglected (this would most likely be the case in central HIC 
where there is no CME). We now assume that the manifest breaking of axial current conservation, 
generated by quark masses, can be neglected for small masses and finite volumes (as there are no 
proper zero modes then). Then the induced axial charge turns out to be conserved provided that 
the topological charge is not changed during Tg re ^ a ^. 

In HIC one can conceive of the following evolution: at the initial stage the nuclear matter 
is compressed and heated (during a time Seating < ^-5 fm) and a topological charge emerges 
due to a sphaleron transition [17]. For light quarks the creation of topological charge leads to 
the generation of an axial charge. At the next stage the nuclear matter cools down (for a time 
Tfireball ~ 5 — 10 fm ) until freeze-out. During this period the topological charge is supposedly 
conserved. The axial charge would be conserved too provided that the quark mass term breaking 
chiral symmetry remained subdominant. The characteristic oscillation time is governed by inverse 
quark masses. Evidently for u, d quarks l/rh g ~ 1/5 MeV -1 ~ 40 fm 3> ^fireball an d the left-right 
quark mixing can be neglected. But it is not the case for strange quarks as l/m s ~ 1/200 MeV -1 
~ 1 fm ^ r hreball anc ^ even ^ a topological charge persists during fireball lifetime, the mean value 
of strange quark axial charge is around zero due to left-right oscillations. 

Thus for u, d quarks, QCD with a topological charge (AT5) 7^ can be equally described at the 
Lagrangian level by topological chemical potential \x$ or by axial chemical potential ^5 

(AT 5 ) * -L(Qg) <=► M 5 c (5) 

simply adding to the QCD Lagrangian A£^ p = /igATs or, alternatively, AC q = /U5Q5. 

In a Lorentz invariant form one can associate a classical background a(x), depending smoothly 
on time, so that 

N f 1 

A£ G = -^K^aix) = -^fj, e K <^=> M5<no75<?; ^5 - a{t) ~ constant. (6) 

Thus, we assume that a time dependent but approximately spatially homogeneous background 
of a pseudoscalar field a{t) is induced at the energy densities reached in HIC during the fireball 
evolution and its gradient defines a 4-vector £ M ~ d^a ~ ((",0,0,0). The field a(t) could be either 
SU(3)f singlet or diagonal part of SU(3)f octet or even a mixture of the two. 

If we want to investigate LPB in HIC with the help of electromagnetic probes, we have to 
account for the photon contribution to the singlet axial anomaly 

Qg -> Q5 = Ql - T 5 em , T| m = ^j^d 3 x e m Tr (i^ fc i') . (7) 
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Now fi 5 is conjugated to the (nearly) conserved Q5 but not to Q5 itself [15] . 

After bosonization of QCD in the light meson sector we use the quark-meson description of 
with vector meson and photon fields appearing in the quark covariant derivative. The anomalous 
PCAC relation ([3]) is bosonized following the Vector Dominance Model (VDM) prescription [191120]. 

The extra term in the Lagrangian takes the form [2T] AC ~ e^f^Tr C,^V v V pa , with = £<5 M o 
for spatially homogeneous and isotropic fireball (the hat denotes the isospin content in matrix 
notation). Therefore a non-trivial topological chemical potential pe is transmitted via axial chemical 
potential to a non-trivial axial charge inducing the photon and meson Chern-Simons Lagrangian 
(a similar motivation can be seen in [2]). And vice- versa, after measuring the LPB background £ 
one could, in principle, find the value of the mean topological charge in fireballs. 

When nuclear matter is superdense, pseudoscalar condensation in the 1 = 1 channel may occur 
too, as it has been outlined in [3] (and references therein). This situation corresponds to introducing 
a vector chemical potential p. As indicated in the introduction, the appearance of LPB is inferred 
in this case from effective Lagrangian techniques, and a direct lattice confirmation is still missing. 
In spite of this, the conclusions of the analysis in [3] seem rather robust: there is a range of densities 
where a LPB vacuum is energetically favorable. Thus for light quarks in hot and dense nuclear 
matter the matrix structure of the Chern-Simons vector £ in flavor space generically includes not 
only SU(3)f singlet but also neutral components of SU(3)f octet, that corresponds to including 
both a fis and a p term, respectively. We note that the appearance of a condensate in the T3 
direction implies in addition a breaking of the isospin symmetry of the vacuum. 



3 VDM Lagrangian in the presence of P-breaking background 



The appropriate framework to describe electromagnetic interactions of hadrons at low energies is 
the VDM containing the lightest vector mesons p°, u and <f> in the SU (3) flavor sector. Quark-meson 
interactions are described by 



Ant = Ql^V^q; 



/It Is; 



(8) 



^ + jrl — 4l s , g w ~ g p = g ~ 6 < ~ 7.8; I q and I s are the identities in the 



where Q - y t 3 
non-strange and strange sector, respectively; and A3 is the corresponding Gell-Mann matrix. The 
Maxwell and mass terms are 
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, det (m 2 ) = 0, 



where (V^ a ) = (A^, uj p , p° = p p , and m v = m 2 p = 2g 2 p fl ~ m 2 . This matrix reflects the 



VDM relations at the quark level [121 [2U] . Finally, in a pseudoscalar time-dependent background 
the Lagrangian contains a parity-odd Chern-Simons (CS) term 



4 



(»V v (x)V pa (x) 



2 



C ejM Vj d k Vi = - C ejki Vj.a N ab d k Vi )b , (9) 
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which additionally mixes photons and vector mesons due to LPB. With this definition, the relation 
( = N c g 2 fj l5 /8ir 2 holds. Numerically ( ~ I.5//5. 

At this point, one should recall the arguments given in the section 2 about the fact that the 
strange quark flavor approximately decouples from the rest due to intensive left-right oscillations 
as m s S> 7fl re bail- ^ n consequence, even if a LPB condensate appears in the strange sector, its 
average will be essentially zero. In addition the role of <j) meson is negligible as its typical mean 
free path ~ 40 fm makes it less sensitive to medium effects. In this sense, from now on we shall 
consider a general £ to be a linear combination of isosinglet and isotriplet cases. Therefore only the 
flavor matrices I q and A3 are relevant and the <j) meson will not be considered in the mass matrix. 

For isosinglet pseudoscalar background C, = I q 2(^/g 2 , the mixing matrix reads 
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, det(AT) = 0. (10) 

SU(2) f 



This case is expected to be the dominant one in experiments where the temperature in the nuclear 
fireball is much larger than the chemical potential T > as it is the case for the LHC or RHIcH. 

Let us now try to estimate the value of £ using simple considerations. A first naive estimate 
would be to assume simply that £ ~ r _1 so a formation time of the pseudoscalar condensate r ~ 0.5 
fm leads to ( ~ 400 MeV. On the other hand, we can estimate ^5 by assuming that the change in 
the free energy has to be of order f n and the average change in the axial charge (AQ5) ~ 1. Then, 
assuming a linear response, ^5 ~ /„- and ( ~ 150 MeV. We conclude that the natural value for £ 
is in the few hundreds of MeV. That is to say that its effects on low-energy hadronic physics are 
potentially large and let us see how its presence affects the in-medium meson dispersion relations. 

The mass-shell equations for vector mesons read 

K%V Vtb = 0; k"V v>b = 0, = g^(k 2 l-m 2 )-km-ie^C P KN e , (11) 

and select three physical polarizations vectors (e + ,e-,ei) for massive vector fields that couple to 
conserved fermion currents (the polarization projectors are described in |15|). The longitudinal 
polarization is orthogonal to k^ and to The mass of this state remains undistorted while the 
transversal polarizations satisfy 



if£4 = [k 2 I - m 2 ± V(C • k) 2 - ( 2 k 2 iV)4; m 2 v± = k 2 - k 2 ~ m 2 v T Q\k\- (12) 



The spectrum is found after simultaneous diagonalization of Tn 2 \su(2) ~ N with ~ (C, 0,0,0) = 
constant. We notice that in the case of pure isosinglet pseudoscalar background massless photons 
are not distorted when mixed with massive vector mesons. In turn massive vector mesons split 
into three polarizations with masses < rny L < m 2 / _. This splitting unambiguously signifies 

parity breaking as well as violation of Lorentz invariance (due to the time-dependent background) . 
Note that the position of resonance poles for ± polarized mesons is moving with wave vector \k\ and 
therefore they appear as broadened resonances, leading to an enhancement of their contribution to 
dilepton production away from their nominal vacuum resonance position. 



Experiments at FAIR and NICA will explore the region /1 ^> T in the future 
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4 Parameters of the VDM Lagrangian and new processes 



The VDM coefficients in ([9]) and (jlOp are obtained from the anomalous Wess-Zumino action [21] 
and related to the phenomenology of radiative decays of vector mesons [22] . In particular the VDM 
coefficients in (jlOp can be estimated from the experimental decay constants [23J for the processes 
7] — > 77, rj' —7- 77, w — > 777, p° — > 777 after taking into account a strong % — 770 mixing [24J. Only 
the ratio of the decay widths uj — > r/7, p° — > 777 is a little sensitive to the mixing and confirms the 
off-diagonal elements of (jlOp , Armed with the above information and after estimating £ one can 
easily determine the modifications on the vector meson spectrum and their eventual relevance for 
the dilepton excess. 

However, this is not the end of the story. If there is parity breaking many other processes are 
possible. They can be estimated by using the spurion technique; i.e. /J5 as the time component 
of a fictitious axial field. Two new processes are then likely to be most relevant inside the fireball 
thermodynamics: the decays 77, 77' — > ittt that are strictly forbidden in QCD on parity grounds. 

To get a rough estimate of the relevance of these previously forbidden processes we take the 
parity even sector and extend the covariant derivative by including a axial chemical potential, 
D u D u — i{n55o u , •} = D y — 2i/j5#o;/. There is no contribution from dimension two operators, 
but dimension four terms (see e.g. [25] ) lead to 

~-^Ld77Tr(d7rd7r); = 77,77', (13) 

where L is a combination of the Gasser-Leutwyler constants £1,2,3 ~ 10~ 3 . To get a numerical 
estimate we take the average pion momenta to be ~ m^/2. The effective coupling constant affecting 
this operator can be estimated to be ~ 0.4, which is large enough to induce substantial 77 meson 
regeneration in the hot pion gas. A very rough estimate of the partial width for the exotic process 
under discussion gives ~ 100 MeV, when we assume ( ~ 200 MeV, to be compared for 

instance to T p ^ n7T ~ 150 MeV. Clearly if p's are in thermal equilibrium in the pion bath, so will 
the 77. A similar analysis leads to an even larger width for the 77' due to the above parity breaking 
operator. On the other hand, the uj decays are not modified by a parity breaking isosinglet spurion 
(but they would for an isotriplet one). 

However at this point, one should be aware that the previous estimate using the chiral La- 
grangian may not be reliable at all because the numerical value of the parity breaking terms is 
quite large. Then one expects substantial mixing with the scalar partners of 77, 77' (i.e. a and /o) 
with comparable masses. A detailed analysis will be given elsewhere. On the contrary we do not 
expect substantial mixing between p and a\ for an isosinglet condensate. 

We have to retain the very important point that a number of light hadronic states will be in 
thermal equilibrium (so their respective abundances will be governed by the Boltzmann distribu- 
tion) and regeneration of these resonances will take place in the fireball similarly to the one taking 
place for /?'s. They will be much more abundant in the fireball if there is LPB and the decays of 
these resonances will be an important source of dileptons so far unaccounted for. 

5 Dilepton production rate in P-breaking medium 

While during the first stage of the HIC, nuclear matter is being compressed and heated (during 
r heating) an< ^ a P ar ity breaking condensate a(t) appears and grows, during the expansion of the 
fireball this condensate dilutes during a time Tfirebajp One can approximate both regimes by a 
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linear function with slopes ~ Cheating an< ^ ~ ^fireball' w ^ n opposite signs. Of course, the cooling 
time greatly exceeds the time taken for the LPB phase to form so most of the expected effects 
should come from this latter period. For the very LPB effect only |£| is important because a 
change of sign causes an interchange in polarizations but does not affect the splitting of masses for 
different polarizations. Naturally there is some error in approximating the time dependence of the 
condensate with a linear parametrization (which is the one where our results are strictly valid) but 
qualitative aspects should be well captured assuming an average or effective constant value for (. 

Dileptons are produced in a variety of decays. The dominant processes in the range of invariant 
masses 200 MeV < M < 1200 MeV are assumed to be the p and ui decays into lepton pairs, the 
u) — > tt°£ + £~ Dalitz decay and similar Dalitz decays for 7] and rj' . These are the basic ingredients of 
the so called 'hadronic cocktail' |10j . conventionally used to predict (unsuccessfully for central HIC) 
the dilepton yielcfl Naturally, the dispersion relations of p, ui, etc. are modified due to conventional 
in- medium effects but these modifications are insufficient to explain the abnormal dilepton yield. 

LPB modifies the calculation in two ways. On the one hand generates new in-medium effects 
on resonances that can be unambiguously predicted in term of the parity breaking parameter £ 
alone. In addition, some hadronic states mix as a consequence of the parity breaking effect and 
the new allowed interactions are able to thermalize some of the lowest lying states, whereas only 
p's could reasonably be expected to be in thermal equilibrium with the pion gas without LPB. In 
what follows we will consider the modifications in the dilepton spectrum due to the p and ui and 
shall postpone discussing the effects of the new processes allowed by LPB. 

The production rate of lepton pairs takes a form similar to the one given in [10] but with 
modified form factors due to LPB, according to our previous discussion 

dN a 2 / n 2 v ml\ 3/2 ^ 1 



d 4 xd 4 kd 2 p T 48vr 2 M 2 V M 2 J \E k p\\ - k\\E 



1 , s m vJ 1 + 



r 2 

1 V 



e 1 ( M 2 -m 2 \ +m 4 ^ 



where ny = 2,0 for p and cj cases respectively, and M > nym n . My is the transverse mass 
My = M 2 + kj, while kx and k\\ are the perpendicular and parallel components, respectively. The 
projectors P<f v are detailed in [15]. A simple thermal average with the Boltzmann distribution has 
been included |26j . T being an effective temperature |9j. Finally, the constants cy normalize the 
contribution of the respective resonances. 

Being a theoretical paper, we will not pursue a detailed comparison with experiments. However, 
in order to make our plots meaningful we have implemented some values relevant for the PHENIX 
experimental set-up. For instance since the oj is not in thermal equilibrium inside the fireball, 
its relative weight c u is normalized to the peripheral HIC result (agreeing with existing pp and 
p- nucleus data 0[9]) and a relative enhancement of 1.8 is assumed for the c p /c w ratio based on the 
data ([9]; V.L. Rykov, private communication). Likewise the integration of (|14p is performed over 
the acceptance region of the PHENIX experiment (pr > 200 MeV and \y e e\ <0.35) and T ~ 220 
MeV. 

Then one is lead to the dilepton production result that is shown in Figure Q] (left), where the 
p spectral function is presented together with the separate contributions for each polarization for 



4 In addition there is a substantial cc background that is not affected by the present considerations. 
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Figure 1: Left: The polarization splitting of the p contribution to dilepton production is shown 
for LPB with £ = 400 MeV. The comparison with the case £ = (shaded region) is presented. 
Right: the in-medium contribution in the p and u channels (solid and dashed line, respectively) 
is presented for Q = 400 MeV together with their vacuum contributions (light and dark shaded 
regions, respectively). The in-medium p yield is enhanced by a factor 1.8 (see text). The vertical 
units are taken to coincide with PHENIX experimental data [9], as well as experimental detector 
cuts and temperature. 



C = 400 MeV and a comparison with Q = (no LPB). Similar results may be obtained modifying 
this parameter, so when £ increases, the circularly polarized resonances appear to be more separated 
from the vacuum one. Therefore a measure of dilepton polarization event-by event may reveal in 
an unambiguous way the existence of LPB, confirming the hypothesis of pseudoscalar condensate 
formation in HIC. 

As already mentioned for central HIC, particularly at low pt, the p/oj production ratio needs to 
be enhanced by a factor 1.8 in PHENIX [9j. This enhancement reflects the multiple regeneration of 
p mesons through tttt fusion in hot pion gas. There is no such a regeneration for narrow resonances 
u and (p. The simulation of u meson production [27] shows that a significant fraction of them decay 
inside of the nuclear fireball and therefore LPB distorted uj mesons may also be responsible for a 
fraction of the abnormal production of dileptons. This is shown in the right panel of Fig. [TJ In 
this plot, the separate contributions of p and u are presented for Q = 400 MeV to be compared 
with their vacuum results. The effect shows a clear increase in the dilepton production away from 
the p — co peak due to the mass shifting of the circularly polarized resonances in the LPB phase. 

Other processes in the region 200 MeV < M < 1200 MeV, especially relevant for dilepton 
production below the p — oj resonance, are the u, r\ and rj Dalitz decays, described by the Kroll- 
Wada formula [28] that includes the contribution of vector mesons and it remains valid in the 
case of LPB provided that we replace the vector meson masses by the values in (|12p according 
to the intermediate meson polarization (L,±). We have checked that this contribution shows a 
substantial enhancement but this and other hadronic processes relevant for dilepton production 
will be discussed in a separate publication. The fact that the decaying meson is distorted by 
the medium complicates analytically and numerically the calculation due to the lack of Lorentz 
invariance generated by a polarization and momentum dependent meson mass. 
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6 Conclusions 



We have explored the consequences of assuming LPB via an isosinglet condensate in HIC. We would 
like to emphasize the simplicity of the approach presented here. The fits presented use the values 
(effective temperatures, normalizations, etc.) quoted by the experiments themselves. The only free 
parameter is £, which is expected to depend on the characteristics of the HIC. It should also be 
said clearly that the presence of LPB does not preclude other many body or in-medium corrections 
[TO! [TH [131 [2], as l° n g as they do not represent double counting. 

A clear signal of LPB would verify that dileptons produced for values of the invariant mass 
above and below the p — oj pole are predominantly of opposite circular polarizations in event-by- 
event measurements. This requires searching for asymmetries among longitudinal and transverse 
polarization for different values of M in event-by-event measurements. A more indirect verification 
would be a detailed account of the dilepton enhancement in HIC. In this paper we have worked out 
in detail the modification of the p and oj dispersion relations, leaving other effects for a forthcoming 
article. 

We have seen that in presence of an isosinglet time-dependent pseudoscalar background the 
vector meson, p and u, propagators are severely distorted. We have computed this effect and found 
that it naturally tends to produce an overabundance of dilepton pairs in the p — oj resonance re- 
gion. We have also shown how LPB induces thermalization in the pion gas of other light resonances 
and how this could enhance the dilepton production. Thus LPB seems relevant to explaining the 
PHENIX/CERES/NA60 'anomaly'. More work is needed before definite conclusions can be drawn 
but we believe that sufficient evidence is accumulated at present to bring these tentative conclusions 
to the attention of the interested readers. 

We acknowledge the financial support from projects FPA2010-20807, 2009SGR502, CPAN (Con- 
solider CSD2007-00042). A. & V. Andrianov are supported also by Grant RFBR 10-02-00881-a 
and by SPbSU grant 11.0.64.2010. 
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